Nonequilibrium critical dynamics of a three species monomer-monomer model 
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We study a three species monomer-monomer catalytic surface reaction model with a reactive 
steady state bordered by three equivalent unreactive phases where the surface is saturated with one 
species. The transition from the reactive to a saturated phase shows directed percolation critical 
behavior. Each pair of these reactive-saturated phase boundaries join at a bicritical point where the 
universal behavior is in the even branching annihilating random walk class. We find the crossover 
exponent from bicritical to critical behavior and a new exponent associated with the bicritical 
interface dynamics. 
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Nonequilibrium models with many degrees of freedom 
whose dynamics violate detailed balance arise in studies 
of biological populations, heterogeneous catalysis, fluid 
turbulence, and elsewhere. The macroscopic behavior of 
these models can be much richer than that of systems 
in thermal equilibrium, showing organized macroscopic 
spatial and temporal structures like pulses or waves, and 
even spatiotemporal chaos. Even the steady state be- 
havior can be far more complicated, involving for ex- 
ample generic scale invariance. Like their equilibrium 
cousins, nonequilibrium systems at continuous transi- 
tions between steady states show universal behavior that 
is insensitive to microscopic details and depends only on 
properties such as symmetries and conservation laws. 

One of the most common continuous phase transitions 
in nonequilibrium models is a transition to an absorbing, 
noiseless steady state jjj, the term absorbing indicating 
the state cannot be left once it is reached. Examples 
of this include directed percolation (DP) the con- 

tact process J4j, the dimer poisoning transition in the 
ZGB model || for the catalytic oxidation of CO, auto- 
catalytic reaction models || , and branching annihilating 
random walks with odd numbers of offspring |?],|| . Both 
renormalization group calculations (|J^] and Monte Carlo 
simulations |j|-|8| show that these models form a single 
universality class for a purely nonequilibrium model with 
no internal symmetry in the order parameter. 

Recently, a number of models with continuous adsorb- 
ing transitions in a universality class distinct from di- 
rected percolation have been studied. These models in- 
clude the probabilistic cellular automata models studied 
of Grassberger et al. |l(]], certain kinetic Ising models 
Pi , the interacting monomer-dimer model [^2|,^3), and 
branching annihilating random walks with an even num- 
ber of offspring (BAWe) RQ. All of these models ex- 
cept for the BAWe have two equivalent absorbing states 
indicating the importance of symmetry of the adsorbing 
state to the universality class. However, the universal 
behavior of this new class is apparently controlled by a 



dynamical conservation law. If fl0|| the important dy- 
namical variables in this class are defects represented by 
the walkers in the BAWe model and the walls between dif- 
ferent saturated domains in the other models, the models 
have "defect parity" conservation law where the number 
of defects is conserved modulo 2. Recent field theoretic 
work confirms this viewpoint p| . 

We study here a simple catalytic surface reaction 
model with three different equivalent monomer species. 
This model could represent either a system with three 
different chemical species or an auto-catalytic reaction 
system in which one chemical species can adsorb on three 
different types of surface sites. This model has adsorb- 
ing transitions to both one and two equivalent noiseless 
states, and therefore is a good model to study the role 
of symmetry in adsorbing phase transitions. The model 
has two fundamental processes: (a) monomer adsorption 
onto sites of a substrate, and (b) the annihilation reaction 
of two dissimilar monomers adsorbed on nearest-neighbor 
sites of the substrate. Here we consider the model only in 
the adsorption controlled limit where process (b) occurs 
instantaneously. We present here results only for the one 
dimensional version of the model. 

Calling the monomer species A, B and C, the param- 
eters in the model are then the relative adsorption rates 
of the different monomer species pa, Pb, and pc, such 
that pa + Pb + Pc — 1- Using static Monte Carlo sim- 
ulations to find the steady state, and dynamical Monte 
Carlo studies described below, we find the ternary phase 
diagram for the model is shown in Fig. [j], where the hor- 
izontal axis corresponds to the relative adsorption rate 
of A and B monomers pab = Pa/(pa +Pb)- There are 
noiseless phases where one monomer species saturates the 
chain occupying the corners of the phase diagram and a 
reactive steady state in the center. There are continuous 
phase transitions from the reactive phase to the satu- 
rated phases. The monomer densities undergo discon- 
tinuous, first-order, transitions from one saturated state 
to another. The points where the reactive phase and two 
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saturated phases meet are bicritical points |16| ] where two 
lines of continuous transitions meet a line of first order 
transitions. We have also constructed a mean field the- 
ory of the model, to be presented elsewhere JT^] , following 
the methods of Dickman pl| . One unusual feature of the 
mean field phase diagram is that the bicritical points lie 
on the edge of the phase diagram if the correlations are 
correct up to single sites or even nearest neighbor pairs. 
Only when the correlations up to triplets of adjacent sites 
are included does the bicritical point appear inside the 
phase diagram, indicating the importance of reproduc- 
ing the correlations induced by large domains of a single 
saturated phase. 



of the exponents can be made by examining the local 
slopes of the curves on a log-log plot and extrapolating 
to infinite times. Then, the effective exponent S(t) is 

-S(t) = {bx{P(tyP(t/b)]/bxb}, (2) 

77(f) and z(t) being defined analogously. Plotting the lo- 
cal slopes versus t" 1 allows us to determine both the 
exponent and the critical point accurately. At the criti- 
cal point the local slope will extrapolate to a constant as 
t -1 — > with a scaling correction linear in t^ 1 [^(J. Data 
taken away from the critical point will have local slopes 
that curve away from the critical point value as t^ 1 — » 0. 




FIG. 1. Phase diagram showing three saturated phases 
(indicated by the letters), and a reactive phase (the unla- 
beled center region). Solid lines indicate continuous transi- 
tions. Dashed lines indicate first-order transitions. 

We have used dynamical Monte Carlo simulations to 
investigate the universality classes of the continuous tran- 
sitions and bicritical points, the critical dynamics of in- 
terfaces between the two symmetric saturated states at 
the bicritical points, the crossover behavior near the bi- 
critical point, and the first-order lines. We use two forms 
of "epidemic" analysis |||l4|,[l9) following the evolution of 
an initial condition chosen close to the saturated state. 
In the first form (defect dynamics) we use an initial con- 
dition consisting of a single vacancy in a saturated phase. 
The second (interface dynamics) starts from two different 
saturated phases separated by a single vacancy. 

From the simulations we find the probability P(t) the 
system does not fall into the saturated state in t time 
steps, the average number of vacancies per run (n v (t)), 
and the typical size of the defect or interface per sur- 
viving run (R 2 (t)). At a continuous phase transition as 
t — > co it is expected that they obey power law behavior 

P(t)~t~ s , (n v (t))~t\ (R 2 (t))~t z . (1) 

Precise estimates of the location of the critical point and 
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FIG. 2. Effective exponents using Eq. ([]) with 6 = 5 for 
the defect dynamics near the critical point at pab = 0.5 on 
the line where the C poisoned phase meets the reactive phase. 
From top to bottom, the 3 curves in each panel correspond 
to pc = 0.3955, 0.39575, and 0.3960, with the middle curve 
corresponding to the critical point. 

Figure ^ shows the effective exponents of the three dy- 
namic quantities near the phase transition to the C satu- 
rated phase at pab = 0.5. Using 10 5 independent runs of 
up to 10 4 time steps at each parameter value, we found a 
critical C monomer adsorption rate of pc — 0.39575(10), 
and the critical exponents are S — 0.16(1), rj = 0.31(1), 
and z = 1.255(15). These values are consistent with our 
expectation that this transition should be in the DP uni- 
versality class, for which the exponents are S = 0.1596(4), 
r] = 0.3137(10), and z = 1.2660(14) We found simi- 
lar exponents for the adsorbing transition at a number of 
other points along the lines separating the reactive phase 
and the saturated states, indicating the transition be- 
tween the reactive phase and any single saturated phase 
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is always in the DP universality class. 
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FIG. 3. Effective exponents for the defect dynamics near 
the bicritical point where the A and B poisoned phases meet 
the reactive phase as defined in (jij) with 6 = 5. From bottom 
to top, the 3 curves in each panel correspond to pc = 0.121, 
0.122, and 0.123, with the middle line corresponding to the 
bicritical point. 

The same kind of analysis at the bicritical point at 
Pab — 0.5, using an initial condition of a vacancy in 
an A-saturated phase, yields a bicritical point at pc = 
p* c = 0.122(1), and very different exponents, given the 
presence of two-symmetry equivalent saturated phases. 
From 5 x 10 5 runs of up to 10 5 time steps we found 
the local slope data shown in Fig. [|, yielding values of 
6 = 0.29(1), r] = 0.00(1), and z = 1.150(15). These val- 
ues indicate that the bicritical behavior falls in the BAWe 
universality class, for which 6 = 0.285(2), rj = 0.000(1), 
and z = 1.141(2) @. 

For pc < p* c along the A-B coexistence line, a sim- 
ilar analysis yields the sub-critical dynamic exponents 
6 s» 0.5, -q ss —0.5, and z m 1. These exponents also de- 
scribe the two species version of the model which can be 
mapped onto the well known problem of the T — one- 
dimensional kinetic Ising model for which these values 
are known exactly ]22] |. 

To further analyze the importance of competition in 
the growth of two equivalent saturated phases at the bi- 
critical point we also studied the dynamics of an interface 
between those two phases. Starting with a single vacancy 
between the two domains, we used two different methods 
to analyze the behavior of the interface. Since there must 
always be at least one vacancy between two different sat- 



urated phases, in the first method we ignore the survival 
probability P(t) and take (5 = 0. We then measure the 
number of vacancies in the interface (n(t)) cx t n and aver- 
age size of the interface (R 2 (t)) oc t z . From 5 x 10 4 inde- 
pendent runs at the bicritical point, each lasting 10 5 time 
steps, we found the other exponents to be r\ — 0.285(10) 
and z = 1.14(2). This type of interface dynamics has 
been used to study the properties of critical interfaces in 
other models in the BAWe class, where similar results for 
r\ and z were obtained OJljj ]. 
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FIG. 4. Effective exponents, as in Fig. ^, for the second 
type of interface dynamics near the bicritical point where the 
A and B poisoned phases meet the reactive phase. 

In the second type of interface dynamics simulations, 
which has not been studied before, the simulation is 
stopped if the interface between the domains has "col- 
lapsed" back to one vacant site. We introduce a probabil- 
ity of avoiding a collapse P{t) oc t~ 5 and corresponding 
vacancy concentrations (n(t)) cx t v and (R 2 (t)) oc t z . 
Figure ^ shows results from 10 7 independent runs each 
lasting up to 10 5 time steps. We find values of 5' — 
0.73(2), 77' = -0.43(2) and z' = 1.15(2). 

Note the value of the dynamic exponent z or z', which 
measures the size of the active region during surviving 
runs, is the same in both types of interface dynamics 
simulations as that measured for the defect dynamics. 
Furthermore, although the exponents 5 and r\ arc differ- 
ent in the three cases, their sum S + r\ (or 5' + 77'), which 
governs the time evolution of the number of vacancies in 
just the surviving runs, seems to be the same. This in- 
dicates a universal nature of the critical spreading of the 
active region for models with two symmetric adsorbing 
states which is independent of whether defect or inter- 
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face dynamics is being considered. A similar result holds 
for some one-dimensional systems with infinitely many 
adsorbing states [ p3| . 

Assuming this conjecture is true, it should be noted 
that simulations using the first type of interface dynam- 
ics, where (5 = 0, yield no information beyond that ob- 
tainable from simulations employing defect dynamics. 
However, simulations using the second type of interface 
dynamics measure an independent dynamic exponent 5' 
which we expect to be a universal number. It would be 
interesting to measure this exponent for other models in 
the BAWe class. 
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FIG. 5. Location of the critical line as a function of dis- 
tance from the bicritical point. The data falls on a line with a 
slope corresponding to the crossover exponent (f> — 2.1 ± 0.1. 

Finally, we measured the crossover from bicritical to 
critical behavior. Near the bicritical point where the A 
and B poisoned phases meet, the boundary of the reactive 
region is expected to behave as (pab — 0.5) cx (pc —Pc)^^ 
where 4> is the crossover exponent [[l6| . We used the dy- 
namical simulation method to accurately determine the 
location of the DP phase boundary between the reactive 
phase and the A saturated phase near the bicritical point. 
From the log- log plot of pab — 0.5 versus pc ~ P*c shown 
in Fig. | we find <j> = 2.1 ± 0.1. 

We have introduced a convenient model to study the 
role of symmetry in the critical dynamics of adsorbing 
phase transitions. We have shown that the universality 
class of the transition changes from DP to BAWe when 
the symmetry of the adsorbing state is increased from one 
to two equivalent noiseless states. Furthermore, we have 
shown that having a symmetry in the adsorbing states 
introduces a richness into the dynamics that is not possi- 
ble if there is a unique adsorbing state. In particular, the 
critical dynamics of the interfaces between two different 
adsorbing states shows a sensitivity to how the dynam- 
ics is defined, and the survival probability of fluctuations 



in the size of the interface from its smallest value is de- 
scribed by a new universal exponent 5'. However, the 
critical spreading of the reactive region, be it a defect 
in a single phase or a domain wall between phases, ap- 
pears to be insensitive to the choice of initial conditions. 
This appears to result from the fact that large reactive 
regions are insensitive to whether the reactive regions are 
bounded by the same or different saturated phases. 
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